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Abstract 

The non-equilibrium dynamics of the kinetic spherical model with a non-conserved order- 
parameter, quenched to T < T c from a fully disordered initial state, is studied at its upper critical 
dimension d = d* = 4. In the scaling limit where both the waiting time s and the observation time 
t are large and the ratio y = t/s > 1 is fixed, the scaling functions of the two-time autocorrelation 
and autoresponse functions do not contain any logarithmic correction factors and the typical size 
of correlated domains scales for large times as L(t) ~ t 1 / 2 . 
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1 Introduction 



The study of collective non-equilibrium behaviour has received a lot of interest in recent years [TJ. One 
particular aspect is the study of ageing which for example may arise if a many-body system is rapidly 
brought out of some initial equilibrium state by the change of external control parameters such as the 
temperature or an external field. If the change of these variables is such that the equilibrium state of 
the system is either at an equilibrium critical point or else in a coexistence phase with at least two 
equivalent but distinct equilibrium states, then one may observe (i) a slow, non-exponential relaxation, 
(ii) a breaking of time-translation-invariance and (iii) dynamical scaling which are the three defining 
properties of ageing systems, see [21 El SI IS] for reviewsjf] While ageing was first studied in glassy 
systems [5J, it has been realised more recently that its essential properties can also be found in simple 
ferromagnets, which may be easier to analyse and still may provide useful clues for the understanding of 
the more complex glassy systems. In particular, in ferromagnets the form of dynamical scaling is quite 
simple since the linear domain size scales algebraically according to Lit) ~ t x l z . Ageing is conveniently 
studied through the behaviour of two-time observables. For the two-time autocorrelation and (linear) 
autoresponse functions, the most simple kind of dynamical scaling behaviour is 



where r) is the order-parameter at time t and at the location r, h is the conjugate magnetic field, 
a and b are ageing exponents herewith defined and fc,R{y) are scaling functions. These scaling forms 
are expected to hold when both times t, s are large and their ratio y — t/s > 1 is fixed. 

The scaling ( 11. ip is generally valid for systems with a so-called simple ageing behaviour, but there 
are exceptions. In certain cases, the above scaling forms are modified by additional logarithmic factors 
which already manifests themselves in the scaling of the linear domain size, i.e. Lit) ~ it/ hit) 1 / 2 for 
a non-conserved order-parameter, and occurs in many systems where topological defects (e.g. vortices) 
play a role, such as planar magnets, frustrated spin systems, liquid crystals or superconductor arrays, 
see [TJ, [HI El OH ttU E21 EEl [HI US, HEJ [T71 [18] . Recently, for equilibrium critical phenomena, Kenna, 
Johnston and Janke [T9l EH] have constructed a systematic theory for these logarithmic factors, based on 
an analysis of the complex zeroes of the partition function, by which they derive scaling relations between 
the exponents describing eventual logarithmic correction factors to simple scaling. They checked the 
results of their analysis by comparing with the many results available for logarithmic contributions 
in equilibrium scaling (see [2T] for a very recent example). On the other hand, for non-equilibrium 
dynamical scaling much less is known on possible logarithmic factors which makes tests of more general 
ideas difficult. For that reason, we consider in this work the ageing and dynamical scaling of the 
four-dimensional spherical model whose logarithmic corrections factors to equilibrium scaling are well- 
established textbook knowledge [22]. While the scaling of its two-time functions has been analysed in 
great detail, starting with [23J, for either d < 4 or d > 4, to the best of our knowledge, the case d — 4 
has never been explicitly studied. 

In the next section, we define the model and briefly recall those elements of the solution which we 
need for our analysis which is presented in section 3. We conclude in section 4. Technical details of the 
calculation are given in an appendix. 

3 Since each of these three properties can be found alone in situations where one would not speak of ageing, it seems 
reasonable to insist on their simultaneous presence in ageing systems. For example, in any system away from its stationary 
state timc-translation-invariancc will not hold. 
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2 Model and formalism 



The spherical model may be defined in terms of real spin variables S r G K. attached to the sites r of a 
hyper-cubic lattice AcZ^ and which obey the mean spherical constraint 



J2S(t,rY=AT (2.1) 
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where M is the total number of sites of the lattice. The dynamics is assumed to be given by the 
stochastic Langevin equation 

dtS(t, r) = AS(t, r) + i(t)S(t, r) + V (t, r) (2.2) 

where A is the Laplacian with respect to r, $(t) is a Lagrange multiplier chosen such that the spherical 
constraint holds and r)(t, r) is a centred gaussian noise with variance {rj(t, r)rj(t', r')) = 2T5(t — t')8(r — 
r'). In writing this, a choice of units was made such that the corresponding kinetic coefficient is set 
to unity. Throughout, we shall assume a fully disordered initial state (where in particular the order- 
parameter vanishes) described by a gaussian variable and the variance (S(0,r)S(0,r')) = 6(r - r')0 

Since the solution of this model is by now standard, see e.g. [23], we shall merely quote those results 
which we shall need for our analysis of the 4D case. One of the central quantities needed for the analysis 
is the function g(t) := exp(2 f*du%(u)), which, as a consequence of the spherical constraint, satisfies 
the Volterra integral equation 

g(t) = f(t) + 2T fdt'f{t-t')g{t') (2.3) 
Jo 

where the auxiliary function f(t) is given by 

fit) = J e- 2 ^ = (e- 4 %(4t)) d < ~°° (STrt)^ 2 (2.4) 

where B denotes the first Brillouin zone, Iq{v) = \ J^d9e~ ucos9 is a modified Bessel function [27J and 
the dispersion relation is, for a hyper-cubic lattice, u)(q) = X^ =1 (2 — 2 cosgj). In particular, the critical 
temperature T c (d) > for d > 2 and is given by 

1 



2T c (d) 



dt [e- u h{At)) d (2.5) 



from which its numerical value is easily found for any d. For example, T c (4) = 6.45438 . . .. Given 
the function g(t), any observable of interest is readily calculated. In what follows, we shall need the 
single-time correlation function C q (t) given by 

(S q (t)S Q '(t)} = (2n) d 6(q + q')C q (t) (2.6) 
where the (discrete) Fourier transforms are defined by 

/« = £ f^)e-^ ; f(r) = ( f q <** (2.7) 
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4 See [24] for a careful analysis of the mean spherical constraint. In the case of a non-vanishing initial magnetisation, 
the fluctuations in the Lagrange multiplier i(t) must be taken into account [25) . This makes the solution of the model, 
even for d ^ 4, a formidable task. In the interest of a relatively simple presentation, we concentrate on the case of a 
vanishing initial magnetisation only, where the formalism at hand with a non- fluctuating j(t) is sufficient [25] . See |25il26j 
for a detailed discussion of the ageing behaviour when, starting from an ordered state, the quench is made to T = T c . 



One then has 
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g(t) 



f2.8) 



where the initial condition C q (0) = 1 was used. In a similar fashion, one finds for the autocorrelation 
and autoresponse functions [23] 



C(t,s) = (S(t,r)S{s,r)) 
5(S(t,r)) 
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(2.9) 



Finally, we shall also measure the relevant length scale L(t) of dynamical scaling, which is also a measure 
of the linear size of correlated clusters, by considering the normalised second moment of the single-time 
correlator 



L\t) 
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J2reA C (^ r ) " C q (t) 
t + 2Tf*dt> (t-t')g(t') 
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(2.10) 



We remark that the response function, and also the correlators C and the typical length scale L at 
criticality, where the thermal term dominates, merely depend on ratios g(t)/g(t') such that the global 
normalisation of the function g(t) will disappear from these physical observables. 



3 Results 



3.1 Solution of the Volterra equation 

In order to obtain g(t) explicitly, one must solve the Volterra integral equation (12.31) . Since its right- 
hand side has a convolution structure, one considers the Laplace transformation g(p) = C(g)(jp) = 
J^°dte~ pt g(t) and finds 



g(p) 



(3.1) 



1 - 2T/(p) 

From standard Tauberian theorems [2H], the long-time behaviour of g(t) follows from the asymptotics 
of g~(jp) as p — > 0. This in turn requires the leading behaviour of f(p) which may be obtained in a 
well-known fashion, see [281 [291 30J , which we outline in the appendix. Together with the known results 
for 2 < d < 4 and d > 4 [23], we have 
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where a = (8tt) 
and C E = 0.5772 . 



a x = -(87r)- d / 2 |r(l 
. is Euler's constant. 



d/2)\ for d ^ 4, a 2 = J t 
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Figure 1: (a) Behaviour of the function g(t) in d = 4 dimensions as a function of time, above, at and 
below T C (A) = 6.45438 .... In (b) the behaviour g(t) ~ t~ 2 for T = 4 < T c (4) and in (c) the behaviour 
g(t) ~ ln(t) _1 at T = T c (4) is illustrated. 



3.2 Scaling of autocorrelation and autoresponse 

We are now ready to discuss the long-time behaviour of our model where the three cases T > T c , 
T < T c and T = T c have to be distinguished. In figure we illustrate the long-time behaviour of 
g(t) for these three cases, obtained by numerically solving eq. (I2.3p [3Tj . for which g(t) clearly shows 
different behaviour. 



1. T > T c . For quenches into the disordered phase, g(t) is obtained straightforwardly by the method of 

residues _ 

_ /(l/r eq )e^ _ 16tt 2 e*^ 

9{ ] " -2T7'(l/r cq ) " ln(r eq ) - C E - 1 (3 ' 3) 
where the finite relaxation rate is obtained from /(l/r eq ) = 1/2T and reads for T > T c r cq ~ 
r ^ c r ^ n t-t • Since the thermal equilibrium correlation length of the AD spherical model behaves 
as [19j [22] ^eq(^) ~ — T c ) _1 ln(T — T c ) and one expects that the relaxation time r eq ~ £g q , this is 
consistent with known results. The two-time functions relax within the finite time r eq to their time- 
translation-invariant equilibrium values. 



2. T < T c . Now, one has by combining (13.11) and (13. 21) . for p small enough 
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(3.4) 



where the exact equilibrium result M 2 q = 1 — T/T c [22] was used. The presence of terms constant in p 
in (13 .4p signal that it is not enough to concentrate on the asymptotic long-time behaviour of g(t) but 
there are also 'global' contributions to g(t) which must be taken into account. 



Formally inverting the Laplace transform in (13 .4p then leads to the following form for g(t), 



if 2 < d < 4 



g(t) ~ T S(t) + < 



m 

eq eq 



where S(t) is the Dirac delta function. Eq. (I3.5P contains two kinds of terms. First, one has the expected 
'regular' asymptotic form, for t — ► oo, g(t) ~ /(t)M~ 4 and second, there appear 'singular' contributions. 
While the response functions follow from the 'regular' terms alone, the 'singular' contributions are 
important for the correct calculation of the correlation functions, as explained in the appendix^ In 
figure [lb we illustrate the AD long-time behaviour g(t) ~ t~ 2 . We observe that no logarithmic factors 
occur, for d — 4. Therefore, in the scaling limit t, s —>■ oo with y = t/s > 1 fixed, we recover the simple 
scaling behaviour 

C(( ' s) = < ( feTTF) ' R(t ' a ^ik-4(^W (3 - 6) 

which could also have been obtained from the well-known scaling functions for d ^ 4 [23J and performing 
the analytical continuation d — > 4. 

3. T = T c . In this case, we have from (13.21) 

i i f -a^p 1 -^ 2 ; if 2 < d < 4 

g(p) = - w + I -a^p-'iCE + lnp)- 1 ; if d = 4 (3.7) 

c I a 2 X P 1 ; if d > 4 

As shown in the appendix, for d = A this leads to the long-time behaviour g(t) ~ 1/ Int. In figure [lb we 
illustrate that this asymptotic form is indeed compatible with the numerical solution of (12.31) . Taking 
into account also the singular parts of g(t), we find the scaling behaviour 



. . T c 1 /hit 1 ..11 Int 1 . . 

C(M) = 8^7Vln7^I ' ^ = IoWV ln^^iy ( 3 - 8 ) 

The calculations are outlined in the appendix. 

Do these extra logarithmic factors imply that simple scaling is modified in the 4D spherical model, 
at least for T = T c ? Indeed, the answer is negative, since in the scaling limit t, s —>■ oo with y = t/s > 1 
fixed and finite, we can write t = ys, hence ^lnt/lns = A/ln(?/s)/ln,s ~ 1 + |ln?//ln s + . . .. Therefore, 
the scaling functions for d = 4 can be obtained from the known expressions for d ^ 4 [23J by analytic 
continuation d — > 4. The fact that the system is at the upper critical dimension d* = 4 of its equilibrium 
critical point only enters into the additive logarithmic corrections to scaling. In this respect, the AD 
spherical models shows a different behaviour from those systems [9l [101 EH H21 [131 [HI [15l [171 [181 132] 
where logarithmic corrections to non-equilibrium dynamical scaling were seen before. 



3.3 Domain size 

In order to understand better the role of the upper critical dimension d* — 4 for non-equilibrium 
dynamical scaling of the spherical model, we now consider the linear domain size L{t) in order to check 

5 In [53] only the 'regular' asymptotic contributions to g(t) are explicitly given, while the effect of the 'singular' terms in 
(|3.5[) is taken into account by non-asymptotic sum rules, with the same end result for the physically observable two-time 
correlations and responses, if d =/= 4. 
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Figure 2: Long-time behaviour of the linear domain size L(t) for (a) T < T c and (b) T = T c . Recall 
that T c (3) = 3.957 ... and T c (4) = 6.454 . . .. 



for the presence of additional logarithmic factors with respect to the naively expected simple behaviour 
L{t) ~ y/t, valid for T < T c and d ^ 4. In eq. (I2.10p we had already related L{t) to the function g(t) 
which fixes the spherical model dynamics. From this, it is already clear that for T = 0, we simply have 
L 2 (t) = 4dt and since the temperature T is generally thought to be irrelevant for quenched to T < T c 
[2], we would expect that this behaviour of L(t) should remain qualitatively correct for all T < T c (d). 

In order to see that this expectation is indeed borne out, we rewrite eq. (12.101) as the ratio of two 
inverse Laplace transformations, since the two integrals can be seen as Laplacian convolutions 



L 2 (t) = Ad 



C-\p- 2 + 2Tp- 2 g(p))(t) 
C- 1 (p- 1 + 2Tp- 1 g(p))(t) 



On the other hand, for T < T c , we have from H3.4[) that g(p) = (2T C ) 1 (1 
quently, we easily find for all T < T c (d), and for all dimensions d > 2 

L 2 (t) = 4dt (1 + 0(1/*)) 



(3.9) 



T/T c )-\l + o{p)). Conse- 



(3.10) 



This is in agreement with our expectation formulated above and also with more indirect conclusions [23] 
drawn from the scaling of the single-time and two-time correlation functions. In figure [2^ we further 
illustrate this by showing that the ratio L 2 (t)/(td) — > 4 in the long-time limit in both 3D and 4D (in 
the 3D case, we have checked explicitly that convergence occurs for times t > 600). 

In a similar way, we now analyse the critical case T = T c (d). The starting point is again eq. (13.91) 
where we now have to insert eq. (13. 7j) . This leads for t large to (see the appendix) 



L 2 (t) ~t x 



2d (hit 
2d 



l)/(lnt-3/2) 



if 2 < d < 4 
if d = A 
iid>4 



(3.11) 



We see that to leading order as t —* oo, we always have L 2 (t) ~ t and that logarithmic contribution 
at most enter into the additive correction of this scaling behaviour. This shows again the difference 



between the AD spherical model and the systems with time-dependent logarithmic scaling of the form 
L{t) ~ (t/lnt) 1 / 2 studied in the literature [3 El [10l [III [121 [131 [IS [IZl [IH] . In figure [2b, we illustrate 
this for the 3D and AD cases. In the 3D case, a nice convergence towards the amplitude L 2 (t) / (td) = 8/3 
is seen, as expected from (13. lip . In the AD case, we extended the calculations up to t = 10 4 and find 
L 2 (t)/(td) ~ 2.06, not too far from the approximate analytical result (13. lip . 

We conclude that in all cases considered, we have found clear evidence that the typical length scales 
as L[t) ~ t 1 ! 2 and that there no evidence for any logarithmic correction factors present. 

4 Conclusions 

We have considered, clS du CclSG study, the ageing behaviour of the four-dimensional spherical model 
quenched to T < T c from a fully disordered initial state. Since the model is at its upper critical dimen- 
sion, we had expected to find modifications of the usual scaling behaviour by logarithmic correction 
factors which is indeed true when looking at the relaxation time r eq as a function of the temperature 
T, for T > T c , and fully consistent with existing knowledge of these corrections at equilibrium [T9], 20J. 
Surprisingly, when studying the time-dependent scaling, both for quenches into the ordered phase and 
for quenches onto the critical point, our main results eqs. (13 .6[ I3.8[ 13. 1 1[) show standard simple ageing 
and we did not see any evidence for logarithmic corrections appearing in the leading scaling behaviour, 
although additive logarithmic corrections to scaling do occur for quenches to T = T c , see eq. (13.81) . We 
have also found a simple power-law scaling L[t) ~ i 1 / 2 , without additional logarithmic factors, for the 
typical length scale, see eq. (13.111) . Surprisingly, this suggests that the fact of a system being at its up- 
per critical dimension should be considerably less relevant for its non-equilibrium dynamical behaviour 
than it is for its equilibrium critical behaviour. Of course, there may be other reasons for a system to 
develop logarithmic scaling, such as vortices [101 [HI E2l 131 HU HSl [T71 [18] or a roughening transition 
of the interfaces between ordered domains |33j . 

Experimentally, observations of the kind made here could be of relevance for the non-equilibrium 
dynamics of systems at a tncritical point (where d* = 3) and which arise for example in diluted magnets, 
meta-magnets or solutions of long polymers at a O-point. It would be of interest to consider a spherical 
meta-magnet in an external magnetic field and then study its dynamics at the tricritical point at the 
upper critical dimension d* = 3. Indeed, non-equilibrium relaxation at a tncritical point was studied 
long ago in the O(n) model by Janssen and Oerding [32]. At the upper critical dimension d = 3, starting 
form a small initial magnetisation, they find to one-loop order the habitually expected logarithmic 
corrections factors, e.g. in the short-time scaling (M(t)) ~ M (\nt)~ a( - n ^ with a{n) = ^f^p^p and 

similar for the correlation with the initial state C{t) = (M(t)M(O)) ~ t~^l 2 (mt)~ aijl \ However, their 
study only considered single-time observables and therefore cannot illustrate directly the compensation 
of logarithmic factors in two-time quantities which was found in the present work. Since a(n) diverges 
as n — > oo, it appears also possible that the n — > oo limit might have special properties. It remains an 
open problem to what extent the results on the two-time quantities in the spherical model at d — d* 
reported here are generic. 
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Appendix 



We outline some details of the calculation for d = 4 whose results were quoted in the main text [31 J . 
For 2 < d < 4 and d > 4 all results can be taken over from 



First, we analyse the leading behaviour of f(p) for p — > 0. The first step is to decompose the integral 
(we set d — 4 from now on) 

POO 

7(p) = / du e-^ u I (Au) d =: J rcg (p) + J sing (p) (Al) 
J o 

into a regular and a singular part by setting J °° = J" + f°°. The integrals are evaluated in the double 
limit i) > 1, p C 1 and pr/ ~ 0(1) kept fixed. The decomposition of the integral into two parts is 
merely a heuristic device in order to arrive rapidly at the singular terms which will be seen for pr] small 
enough; of course the end result for f(p) should be independent of r]. First, one has, to leading order 
J ve Jp) -> (2T C (4))- 1 . Second, we use the identities eqs. (3.351(4)) and (8.214(1)) from [32] 



duu~ 2 e- pu = pEi(-prj) +T 1 ~ l e~ m (A2) 

Ei(-x) = Cs + lnx + V^TTr- 5 for x > (A3) 

fe=i 

where Ei(x) is the exponential integral and then obtain, to leading order in the double limit introduced 
above 



pEi(-pri) + 



in)- 2 



6 -pv 



V 

~ (8Tr)- 2 [pC E +p\np\ (A4) 

where we dropped an //-dependent term ~ p In rj which should cancel anyway against a corresponding 
correction coming from the regular term f Teg (p). Then the case d — 4 of (13. 2p follows. 

Next, we derive g(p) for d — 4 and T < T c . From the above and A3. If) . the leading terms for p — > 

are 

L_ + C ^ + ^ ln P f A5) 

and g(t) follows by inverting the Laplace transformations, using the identities £ _1 (l)(t) = 5(t) and 

£-\p)(t) = -n / dpe* = *'(f) 

^ J c—ioo 

d 2 r +io ° , e p *lnp d 2 _ 2 



C-\p\np){t) = —J dp—^ = — (InCfe-lnt) = ^ 2 (A6) 

(see eqs. (17.13.1) and (17.13.12) in [34]). Hence we arrive at the expressions ( 13.41) and ( 13. 51) for d = 4. 
From eq. (12.91) we get for the correlation function: 

fd/4 d/4 ps //-(-<? \ 

C(t, .) = <7^y^ + 2TM c 4 q (t,) d / 4 jf dfc / - hj g{t x ). (A7) 

in which we estimate the value of the integral by developing the / function at the first order: 
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2 / 2T c M e 2 q V 2 J t + s M* 
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In the second line we used the sum rule j^°dug(u) = (2T c M eq ) _1 [231 (eq. (2.39)]), which follows in 
our presentation from the first singular term ~ 8(t) in ( 13. and clarifies the important contributions 
coming from the singular terms in g(t). Furthermore it follows from the explicit form of g(t) that the 
regular long-time approximation for g(t) may be used in the second integral in the second line. In 
addition, we replace the derivative of the function / by its asymptotic value. In the last line, we see 
that the second term is negligible in the s — * oo limit. 

Therefore, using the exact expression M 2 = 1 — T/T c , we find for any dimension d 



{t + sfj V eq T c 



eq 1 rp eq 

(t + s) 2 J eq \(y + i) 



where the last expression gives the scaling limit t, s — > oo, with y = t/s fixed. 

The autoresponse function is found directly from (12. 9p by inserting the long-time behaviour for g(t) 
from (I3.5p . 

For the critical case T = T c at d = 4, we have 

167r2 1 1 , A1 „ 

and we now use an approximate method based on the relation lnp = lim ra ^ (p n — l)/ n - 111 order to find 
g (t) , at least approximatively, we write for p small enough 

hm— ^— — -~hm — ^ h . . . (All) 



C E p + plnp n-*o n C E - 1 (1 + p n /(nC E - 1)) n->o nC^ - 1 \j? nC E - 1 
Carrying out the inverse Laplace transform, we find 

£-! f - ^ (t) ~ lim - f 1 — + 

\C E p + plnp) y } ~ n-?anC E -l\ (nC E - 1)T(1 - n) 

n ( e~ nlnt 

~ lim 1 H h . . . 

n^o nC E -l\ (nC E - l)r(l - n) 

1 

In* 



(A12) 



and hence we expect #(t) ~ (hit) 1 for large times. The approximate derivation of (1A12|) relies on com- 
muting several limits. In order to check whether this is justified, we compare in figure [Tb our analytical 



result with the direct numerical solution of the Volterra equation and find that the main feature, namely 
the logarithmic dependence of g(t) on t, is correctly reproduced, but the corresponding amplitude is 
not. In addition, since at criticality only ratios g(t)/g(t') enter into the leading contributions of the 
physical observables of interest, the corresponding amplitudes cancel and hence will not be required. 

Finally, we prove (13. lip . For d > 4, we consider the form (13. 9p together with g~(p) = a 2 1 p~ 1 . Then 

L'jt) = U C ;Y[ + m "T2!f\ 2d t2 =2 dt (A13) 

w C- 1 (p- 1 + (2TJa 2 )p-^(t) (2T c /o 2 )-( ^ ' 

since the second term in both the numerator and the denominator is the dominant one for p —>■ 0. 
Similarly, for 2 < d < 4 we have 

_ C-\p-- + {2T c /a l )p-^){t) ^ l/T(d/2 + l) _ 

L {t) - M c-i(p-i + (2T c / ai )p-^)(t) ~ M i/r(d/2) t ~ 8t (A14) 

Last, but not least, for d — 4, we must find the inverse Laplace transforms of hk(p) '■= p~ k (CE + \np)~ 1 , 
with k — 2, 3. Applying the same procedure we used in ( 1A12I) in order to derive g(t), we obtain 

c-\Hv)W) - ~, . c-\hMM ^ - (M5) 

which completes the proof. We expect that this approximation should give the correct leading t- 
dependence. Since L 2 (t) is given by the ratio of two such expressions, any inaccuracy in the associated 
amplitudes should largely cancel. 



Note added in proof: after this work had been completed, we became aware of the paper which 
studies the ageing in the spherical model with competing interactions. In particular, a detailed analysis 
of g(t) is presented, the results of which are in agreement with our analytical findings, in particular 
g{t) ~ ljhxt at T = T c . 
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